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Abstract 

We use numerical integrations to study the asymptotical 
behaviour of a homogeneous but anisotropic Bianchi type IX 
model in General Relativity with a massive scalar field. As 
it is well known, for a Brans-Dicke theory, the asymptotical 
behaviour of the metric functions is ruled only by the Brans- 
Dicke coupling constant loq with respect to the value —3/2. 
In this paper we examine if such a condition still exists with 
a massive scalar field. We also show that, contrary to what 
occurs for a massless scalar field, the singularity oscillatory 
approach may exist in presence of a massive scalar field having 
a positive energy density. 

Keywords: Anisotropic Bianchi type IX model - Scalar field - Chaos 
- Asymptotical behaviour 



1 Introduction 

Historicahy, scalar fields in cosmology were introduced by Dirac ^ 
and Kaluza-Klein works and then from Jordan, Brans and Dicke 
1^ first scalar-tensor theory. Today they find new justification 
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through particles physics theories. Hence the Higgs mechanism ex- 
plaining the mass of the particles is a massive scalar field [3]. The 
supersymmetry concept which assumes a new symmetry between the 
fermions and the bosons, predicts that several scalar fields should ex- 
ist. In cosmology, the recent acceleration of the Universe expansion 
can be explained by a quintessent scalar field. It is also one of the 
best mechanism to produce inflation at early times. It could even 
mimic the galactic dark matter (HI Ej responsible for the rotation 
curves flattening. 

What about the geometrical framework of this paper? Most of times, 
an isotropic and homogeneous Universe, and thus a Friedmann- 
Lemaitre-Robertson-Walker(FLRW) metric, isl assumed. This is 
nicely confirmed by WMAP observations. However the geometry 
could be different at early times. A way to generalise it is to con- 
sider a homogenous but anisotropic cosmological model of Bianchi. 
Nine Bianchi models exist and the Bianchi type IX one is probably 
the most interesting. Following the Belinskii-Khalatnikov-Lifshitz 
(BKL) conjecture its oscillatory singularity approach would 

be shared by the most general inhomogeneous models 
In this paper, we examine the metric functions asymptotical be- 
haviour of the Bianchi type IX model in General Relativity with a 
massive scalar field. A similar study was performed in .12_ with a 
specific form of potential. With only a perfect fiuid, the metric func- 
tions oscillate in the vicinity of the singularity, adopting the so-called 
Mixmaster behaviour. With a massless scalar field, V. A. Belinskii 
and I. M. Khalatnikov|13j shown that the oscillations were destroyed 
when the Universe collapse to singularity. In fact their results applied 
for a massless scalar field having a positive energy density or, equiv- 
alently, respecting the weak energy condition^. However, this is not 
always the case, in particular if we consider the Brans-Dicke theory. 
Then in the Einstein frame, the behaviour of the metric functions 
strongly depends on the control parameter called the Brans-Dicke 
coupling constant. As it is well known, when lv > —3/2, the singular- 
ity approach is monotonic. The metric functions adopt a kasnerian 

^The weak energy condition is that the sum of the scalar field pressure and 
density be positive. For General Relativity with a massless scalar field for 
which = p^, equation of state of a stiff fluid, it is equivalent to write that the 
scalar field energy density is positive. 
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behaviour, decreasing as some power of the proper time. However, 
when uj < —3/2, the metric functions oscillate forever and there is 
no singularity. Heuristic argument about these facts using the BKL 
approximation are shown in the appendix whereas some considera- 
tions about chaos may be found in or in section El They are 
valid whatever the massless scalar field theory, i.e. even if uj = uj{(j)) 
instead of being a constant. Note also that General Relativity with 
a massless scalar field is equivalent to General Relativity with a stiff 
fluid having a positive or negative energy density. Hence, dynamical 
properties of the Universe (that is singularity, oscillations, monotonic 
behaviour, etc) in such a theory do not depend on the initial condi- 
tions. Moreover, a massless scalar field seems to kill the singularity 
oscillatory approach which is one of the main interest of the Bianchi 
type IX model. 

One of the goal of this paper will be to know what happens in pres- 
ence of a massive scalar field. Considering several scalar field poten- 
tials, we will show that dynamical properties of the Universe now 
depend on the initial conditions. However some properties are still 
ruled only by some control parameters (i.e. the Brans-Dicke cou- 
pling constant and the constant appearing in the potential of the 
scalar field). Using a toy model, a second goal will be to show that a 
massive scalar field having a positive energy density does not always 
prevent from having a singularity oscillatory approach. 
The organisation of the paper is as follows. In section |21 we write the 
field equations. In section |21 we numerically analyse the asymptoti- 
cal behaviours of three theories with respect to their control param- 
eters. In section^ we try to answer the question: is an oscillatory 
approach of the singularity always killed by a massive scalar field 
having a positive energy density? We will produce a toy model that 
shows that the answer to this question is negative. In section El a 
Hamiltonian approach is used to detect the presence of chaos for the 
theories of section O We conclude in section El 

2 Field equations 

The action of General Relativity with a massive scalar field is: 

S = j{R-{oj + 3/2)0-V,/.<A''' - U)^ (1) 
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uj is the Brans-Dicke function describing the couphng between the 
metric and the scalar field (p. U is the potential describing the self 
coupling of the scalar field. 

The diagonal metric of the Bianchi type IX model is given by: 

(2) 

The are the 1-forms defining the Bianchi type IX model. The field 
equations are got by varying the action with respect to the metric 
functions and the scalar field. In the r time such as dt = e°^~^^^^dT, 
they write: 
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A dot means a derivative with respect to r whereas C/^ is the po- 
tential derivative with respect to the scalar field. In this paper, we 
will choose uo = loq where loq is a constant. In the next section, we 
integrate numerically these equations. 

3 Numerical study of the Bianchi type IX 
model 

Without the potential and lo = loq, the above action is the Brans- 
Dicke action in the Einstein frame jl^. Metric functions asymptot- 
ical behaviour (that is oscillation or monotonicity) is then ruled by 
the value of ujq with respect to —3/2 as explained in the introduc- 
tion. However we have to note that when ui < —3/2 there are some 
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known problems. Hence, a scalar field negative kinetic energy leads 
to an instability of the empty Minkowski spacetime or some ghosts 
and tachyons propagation in linearised field equations. 
In this section, we want to know if some similar results exist when we 
consider a massive scalar field: does the metric functions asymptot- 
ical behaviour continue to depend only on some control parameters 
or do they also depend on the initial conditions? 
We will consider three potentials that are widely used to explain 
dark energy or inflation: 

• U = 2A, i.e. a cosmological constant which could model vac- 
uum energy. 

• U = (j)^ , i.e. a power law potential whose Ratr a- Peebles [16] 
one is a subcase with /c < 0. It appears in models of super- 
symmetric QCDpTj. 

• U = e'"^, i.e. an exponential law potential. It can not model 
quintessence but is the outcome of compactification of higher- 
dimensional theories. 

In the Brans-Dicke frame, the two first theories may be cast into the 
Brans-Dicke theory with a power law potential. The third one has 
been studied in ^1^112011^ from the Bianchi models isotropisation 
point of view. In ^21; the case with U = m'^tf)'^ was studied. 
Since we want to find what rules the metric functions behaviour 
(the control parameters? the initial conditions?), we need to identify 
what are these behaviours. During the numerical simulations (see 
below), we identified several of them, sometimes different from those 
encountered in the Brans-Dicke theory: 

• The first one, represented on graph Ai, is such that the os- 
cillations never end. It is what occurs for the Brans-Dicke 
theory in the Einstein frame when the massless scalar field en- 
ergy density is negative (w < —3/2) or equivalently in General 
Relativity with a stiff perfect fluid having a negative energy 
density. Contrary to General Relativity in the empty or with a 
perfect fluid of dust or radiation, these oscillations do not lead 
to a singularity: the Universe 3-volume does not go to zero in 
the r or i time. Moreover, the oscillatory behaviour can not 
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Fi gurc 1; The three metric functions asymptotical behaviours encountered during 
numerical simulations: oscillation(Ai), singularity monotonic approach(Kasnerian for 
-Bi, non kasnerian for B2), expansion (anisotropic expansion for Ci and C'n respectively 
in the t and t times, isotropic expansion for C2 and C21 respectively in the t and t times). 

be continuously approximated by the BKL approximation^ be- 
cause the dynamics has local maxima in the neighbourhood of 
which none of the curvature term may be neglected. 

• The second one is represented on the Bi and B2 graphs. It is 
a monotonic approach of the singularity. 

It may be kasnerian as on Bi graph. It means that the func- 
tions a, (3 and 7 asymptotically behave as some linear decreas- 
ing functions of the r time. This is due to the fact that, in 

^See appendix for a description of this approximation. 
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the isotropy vicinity, the curvature terms and the potential 
become neghgible with respect to the second derivative a, f3 
and 7. The Bianchi type IX model is then approximated by a 
Bianchi type / one. This singularity approach is the same as 
the one of the Brans-Dicke theory in the Einstein frame when 
the scalar field energy density is positive(w > —3/2, see ap- 
pendix). 

It may also be non Kasnerian (but always monotonic) as on 
the figure B2 where the singularity is reached for a finite value 
To of the T time. This singularity also exists in the proper time 
t for a finite value t = to. Once again, the curvature terms may 
be neglected and the Bianchi type IX model approximated by 
the Bianchi type / one. However the potential is not negligi- 
ble any more. It implies that a, /3 and 7 are not some linear 
functions of r. Thus the solution (which in general can not be 
calculated analytically) is not kasnerian. In particular this is 
the case with the power law potential U = (j)'' with k a negative 
integer and </> < (see tabled for a numerical application). 

• The third one is shown on graphs Ci and C2. The oscillations 
end for a finite value tq of r for which the Universe expands 
and the metric functions first derivatives diverge. The met- 
ric functions behaviours of these two figures are not physically 
equivalent. Indeed, in the proper time t, the graph Ci trans- 
forms into graph Cn. We recover a singularity since the Ricci 
scalar diverges. However the graph C2 casts into the graph C21 
where the singularity disappears. In this last case, the slopes 
of the three metric functions tend to the same value and the 
Universe isotropises. 

Note that we have not detected any cigar or pancake singularity. 
It does not mean that they would not exist for other theories with 
some massive scalar fields. Moreover, whatever the functions to and 
U, there is no isotropic point singularity. It is because, the difference 
between two first derivatives of a, /3 or 7 when the curvature terms 
are neglected is always a constant. It is why it does not exist for 
the B2 figure in the r time a similar behaviour as C21 in the t time. 
Some numerical examples of the seven metric functions behaviours 
illustrated on figure ^ are given on table ^ 
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Table 1: Initial conditions of some numerical applications illustrating the 7 metric 
functions behaviours of the figure^ ET means that integration is to early times whereas 
LT means to late time. 



Let US explain how we numerically solved the field equations to get 
the figures I3I6I where the above behaviours were observed. Each 
theory we consider is defined by a couple of control parameters ap- 
pearing in the functions uj, here taken as a constant uq, and U. For 
a constant potential, it is (wq, A) and for a power or exponential law 
potential it is (ujQjk). Each couple defines a plane. We integrated 
and determined the asymptotical behaviour of the metric functions 
for each of their points at early and late times. Numerical integra- 
tions of the field equations have been performed using a Runge-Kutta 
order 5 algorithm. 

What about the initial conditions? If we look at the equation (j7|l. 
we see that it can be rewritten as: 



i+^<^ + e2("+/3+7)C/^ 



3 + 2u; 



Hence we have used for each graph of figures IHIfil two sets of ini- 
tial conditions (ao) "Oi /?0) /^Oi 70i 70i (^o) corresponding to the ± sign 
above. Moreover since each numerical integration uses different val- 
ues of the control parameters, it starts with a different value of (pQ. 
To be sure that the properties observed on the figures are stable, 
we have used several sets of initial conditions. 

Now, let us explain the meaning of the graphs on the figures I3I6I 
Basically each of them represents a plane with ujq as ordinate and 
a second control parameter (A or k) as abscissa. On each figure, 
the left (right) graphs represent the numerical integration results to 
early(late) times. The meaning of the textures covering the graphs 
is the following: 

• The horizontal lines (first texture on figure ^ will represent 
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the points of the plane where a monotonia approach of the 
singularity arises. It may be kasnerian as on graph Bi or not 
as on graph B2 and the energy density is always positive. We 
have not distinguished these two cases because a numerical test 
will be too uncertain to discriminate them. 

• The vertical lines (second texture on figure|2) will represent the 
points of the plane where an oscillatory behaviour of the metric 
functions occurs with a negative energy density. It correspond 
to graph Ai on figure ^ 

• The tilted lines (third texture on figure ^ will represent the 
points of the plane corresponding to an expanding Universe. It 
occurs for a finite value tq of the r time, singular or not, with 
a positive energy density. It corresponds to graphs Ci and C2 
on figure n 

• The points (fourth texture on figureEj) will stand for the points 
of the plane where the asymptotical behaviour has not been 
clearly determined. Normally they should be horizontal line. 
However the monotonic approach of the singularity generally 
arises too slowly and an underflow occurs. 

• The white area will be the points of the plane for which the 
numerical tests have not been able to determine in an unques- 
tionable manners, the asymptotical behaviour before an over 
or underflow. It is not possible to say anything about them. 




Figure 2: The four textures used in the figures r{l6l to represent the asymptotical be- 
haviour of the Universe: contraction with a positive energy density, oscillatory behaviour 
with a negative energy density, expansion with a positive energy density, probably con- 
traction with a positive energy density. 

In the next three subsections, we describe the results of the numerical 
integrations for each potential. 
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3.1 Cosmological constant U = 2h. 

Results of the numerical integrations are shown on figure El The 
abscissa stands for A values and the ordinate for cj values. The metric 




Figure 3: Asymptotical behaviours of the Bianchi type IX model when u) = u)o (in 
ordinate) and U = 2 A (A in abscissa). The left graphs correspond to early times and 
the right ones, to late times. 

functions asymptotical behaviours partly (but weakly) depend on a 
A = Aq critical value related to the initial conditions. Summarising 
these behaviours, we get: 

• When u; < —3/2 and A < 0, the metric functions oscillate 
forever without tending to a singularity. 

• When uj < —3/2 and A > 0, the Universe is spatially open: 
the metric functions are always asymptotically expanding for 
a finite value tq of r and have a minimum. 

• When uj > —3/2 and A < 0, the Universe is spatially closed 
and asymptotically tends monotonically to a singularity, kas- 
nerian or not, at early and late time. Metric functions have a 
maximum. 

• When UJ > —3/2 and A > 0, it always exist a singularity at 
late(early) times. Then, at early (late) times, if < A < Aq, 
the Universe also tends monotonically to a singularity and is 
thus spatially closed. If Aq < A, the Universe is expanding. 
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It is possible to partly explain this distribution of the asymptotical 
behaviours with respect to the control parameters ujq and A. From 
the Klein-Gordon equation, it is easy to show that (f) behaves as an 
exponential of r. Thus (/)^(/)~^ is a positive constant. Consequently, 
when curvature terms are negligible in the constraint equation 
the term with the Brans-Dicke constant dominates. Then as for the 
Brans-Dicke theory, the metric oscillates when uj < —3/2 because 
aP + 07 + /37 < 0. On the contrary when uj > —3/2, the three first 
derivatives may have the same sign and thus the metric functions ap- 
proach to singularity is monotonic. When the Universe expands, the 
numerical simulations show that the dominant terms in the equation 
© are the first derivatives of the metric functions and the potential. 
Consequently, when the three metric functions expand, all the first 
derivatives have the same sign and A must be positive. This is what 
occurs on figure [31 

The asymptotical behaviour of the metric functions is thus weakly 
related to the initial conditions but to the signs of 3 -|- 2u; and A. 
This is due to the fact that the potential does not depend on the 
scalar field. In what follows, we are going to consider some poten- 
tials depending on (/>. Then we will check that the metric functions 
asymptotical behaviour strongly depends on the initial conditions 
although some properties, such as oscillations, stay partly ruled by 
the control parameters. 

3.2 Power law of the scalar field U = cj)'^ 

Numerical results are shown on figure |1] with k in abscissa and lo 
in ordinate. We begin to consider non integer values of k. Conse- 
quently, we choose a positive initial value for the scalar field. This 
guarantees that the potential is real. 

One more time, when uj > —3/2, the metric functions never oscil- 
lates. They monotonically decrease to a singularity or increase at 
late or early times depending on the initial conditions. Some solu- 
tions have both a Big-Bang and a Big-Crunch and are thus spatially 
closed. Others begin(end) with a singularity and end(begin) by an 
expansion. In any case, the singularity approach is monotonic. Solu- 
tions do not exist that are spatially open at both early and late times. 
The set of parameters couples {uj, k) corresponding to an expanding 
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Universe has a hyperbolic form and is centred on the ordinate axis. 
It is explained in part by calculating the necessary condition such 
that the Universe becomes isotropic and expanding. Using [22j, it 
writes lo < ^(/c^ — 3) and is valid when 3+2w and U are positive. It is 
the equation of an hyperbole in the (w, k) plane with a minimum in 
{uj,k) = (—3/2,0) in agreement with the figure IH Note that all the 
points of the hyperbole do not necessarily correspond to an isotropic 
Universe as explained in section |31 with Ci and C2 graphs. 
When Lo < —3/2, the metric functions oscillate or increase at early 
or late times, depending on the initial conditions. Some solutions 
have a minimal 3- volume. There is no solution which asymptotically 
oscillates at early and late times contrary to what happens in pres- 
ence of a cosmological constant or for a massless Brans-Dicke theory. 
It shows that the positive potential ends to dominate the negative 
kinetic term of the scalar field. It thus destroys the oscillations and 
leads to a monotonic behaviour. 

What happens when k is an integer? The asymptotical behaviour 
of the metric functions then partly depends on the k parity as shown 
on figure El In order to find some results different from the previous 
one, we choose a negative initial value for the scalar field. The two 
upper graphs are such that the root of (j)'^ is positive, and the two 
other ones with the negative root. 

When u)Q > —3/2, again, there is no oscillation. The Universe is 
expanding or monotonically approaches a singularity. This is not an 
obvious result, because the negative potential could dominate the 
scalar field kinetic terms thus implying some oscillations. However, 
it does not occur. We see that an expanding Universe only appears 
for even values of k corresponding to a positive potential as above 
when we considered real values of this parameter. As opposed to this 
case, the monotonic approach (kasnerian or not) of the singularity 
does not depend on k being even or not. 

When ujQ < —3/2, the Universe is either expanding or oscillating. 
It never contracts to a singularity. Once again, expanding Universe 
only appears for even value of k, i.e. positive potential, whereas 
oscillating behaviour of the metric functions does not depend on k 
parity. However, stable oscillating solutions at both early and late 
times only appear for odd values of k whereas some solutions having 
a minimum 3-volume only arise for even values of this parameter. 
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This application shows in an obvious way the dependence of the Uni- 
verse asymptotical behaviour on the sign of the scalar field and thus 
of the potential. 

3.3 Exponential law of the scalar field U = e'"^ 

The results of the numerical integrations are shown on figure IHl with 
k in abscissa and in ordinate. These figures have been computed 
assuming (pQ < 0. Apparently the metric functions oscillations only 
occur for < at early times. But if we choose a positive cpo, it 
also arises for A; > and again at early times. We have not observed 
such behaviours for late time. Only two properties seem indepen- 
dent from the initial conditions. The fist one is a singularity at early 
times when to > —3/2. It is reached monotonically by the metric 
functions. The second one is the Universe expansion at late times 
when CO > —3/2 and k > 0. 

For the other sets of uj and k values, the asymptotical state depends 
on the initial conditions. 

When u! > —3/2 and /c < 0, it exists some solutions having a Big- 
Bang and Big-Crunch whose approach is monotonic. No oscillatory 
behaviour has been detected. Some solutions with a Big-Bang at 
early times and which expands at late time have been detected but 
never the opposite. There is no solution which is expanding both at 
early and late times and thus has a minimum 3-volume. 
When UJ < —3/2 and k < 0, the early and late times are such that 
the Universe is expanding or oscillating. Whatever k, it exists some 
solutions having a minimal 3-volume. It does not exist any solution 
oscillating at both early and at late times, indicating that this state 
is unstable. 

To conclude, this section shows that, contrary to what happens with 
the Brans-Dicke theory in the Einstein frame, the asymptotical be- 
haviour of the metric functions also depends on the initial conditions 
and not only on the control parameters. However, it is always pos- 
sible to find some properties which are independent from the initial 
conditions, in particular with a cosmological constant. As instance, 
forever oscillating Universe or Universe having a minimum volume. 
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4 Positive energy density and metric oscilla- 
tions 



One of the interest of the Bianchi type IX model is its singularity 
oscillatory approach |23j. It is able to isotropise the Universe or iron 
out inhomogeneities. Such a singularity exists in General Relativity 
in presence of a perfect fluid of dust or radiation whose effects are 
then negligible. In presence of a massless scalar field, things are dif- 
ferent. When the energy density is positive or equivalently the weak 
energy condition is respected {u > —3/2), the singularity approach 
is kasnerian and monotonic. Otherwise(u; < —3/2), the metric func- 
tions oscillate as on graph Ai but the Universe does not go to a 
singularity (see appendix for heuristic arguments with the BKL ap- 
proximation or the next sections). Apparently, the situation could 
be the same with a massive scalar field: the three basic theories we 
worked out in the previous section always show a monotonous (kas- 
nerian or not) approach of the singularity and never an oscillatory 
one. If it was true, it would mean that a scalar field, massless or 
massive, kills the oscillatory approach of the Bianchi type IX sin- 
gularity. Here, we show using a toy model that this assumption is 
wrong: 1 - singularity oscillatory approach may exist in presence of 
a massive scalar field. Moreover, 2 - it does not require a negative 
scalar field energy density but the non respect of the weak energy 
condition could he necessary. This difference between "a negative 
scalar field energy density" and "the non respect of the weak energy 
condition" only exists for a massive scalar field. For a massless one, 
there is no difference. The potential thus plays an important role. 
To show these two important results, first we write the density and 
pressure of the scalar field as: 





and we define the following variables: 
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K2 = y{pcp-Pcp)p^^ 

V being the comoving 3- volume e""'"^"^'^. Ki has no particular mean- 
ing but is a positive quantity. When P4,/ tends to a constant, K2 is 
proportional to the variation of the comoving 3-volume in the proper 
time. Using Ki and K2, the field equations write: 
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ki = K2K1 



We then define our toy model as Ki = kie^ > Q where ki is a 
positive constant. The scalar field energy density is thus positive and 
corresponds to K2 = 1, i.e. the comoving 3-volume asymptotically 
behaves as t. We numerically integrate the field equations to get the 
two graphs of the figure [Tj The first one is obtained by integrating 
to early times. Then, the scalar field energy density tends to vanish 
and the metric approach a singularity by oscillating. The second 
one is got by integrating to late times, and the comoving 3- 
volume increase as the metric functions oscillate faster. It allows to 
the metric functions derivatives to counterbalance the increasing of 
/9(^y2 which stays negligible in the field equations. It is not possible 
to recover analytically the forms of u) and U . However, we can 
check numerically that uj is not a constant and is always smaller 
than —3/2. It means that it is the potential which keeps the energy 
density positive all along the integration and that the weak energy 
condition is violated. 

Hence our toy model allows an oscillatory singularity approach with 
a massive scalar field having a positive energy density but violating 
the weak energy condition. 
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5 Hamiltonian analysis 



In this section we would like to use a simple method to foretell the 
presence of chaos or not in the various systems that we studied above, 
in particular by making predictions in the {ujq, k) plane. 

5.1 Chaos energy threshold 

We shall use an approach of chaos in term of energy threshold. It has 
been shown by several methods (numerically and analytically, see for 
example a recent study by ) that such a threshold of chaos exists 
in the Bianchi type IX vacuum Universe case but also in presence 
of massless scalar fields. 

In the vacuum case, the chaos threshold is reached when the Hamil- 
tonian density 

/i(a,/3,7) = Cc + Cp 

with 

Cc = d/3 + 07 + /37 

+ + e^^ - 2(e2("+'3) + e2("+^) + 6^^"+"'^) 

is vanishing. This threshold is precisely realized since one has ex- 
actly h = 0. 

In presence of ordinary matter, the curvature terms dominates the 
matter ones that behave as pV'^. Asymptotically there is a singular- 
ity implying a vanishing volume V. Thus we get again in this limit 
h = with a chaotic approach to the singularity. 
In the case of stiff matter, the matter component adds a contribu- 
tion pV'^ that scales like poV^'^~^') with 7 = 2. Thus the constraint 
equation at any time is now h — po = 0. For /jq > we find that the 
chaos is stopped since h > (we stay below the threshold). But for 
an exotic matter such that po < we get h < and we enter in a 
chaotic regime. 

The same conclusion applies in the case of a massless scalar field (f), 
equivalent to a stiff matter. It adds a kinetic contribution to the 
Hamiltonian density such that the constraint reads now h=^^(w + 
|). Hence for 1^ > — | there is no chaos but there is for u < — |. 
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Other cases have been studied including both ordinary matter and 
a scalar field, leading to the same conclusion jl4j. 



The problem here is to predict what is going on with a massive scalar 
field, but without ordinary matter. The corresponding Hamiltonian 
density for the scalar field reads : 

Hi4^) = E,{4>)+Ep{(t>) 

with 

EM = -\VV^ = -ie=<"+''+^lf/ 

and the constraint for the total system is thus 

h + H{<i)) = (8) 

From equation © chaotic situations can be predicted in several cases 
since here the chaos condition reads h=-H((/)) < 0. Simple predic- 
tions will be made in presence or in absence of volume singularity 
(see next section) and are summarized in the section entitled general 
predictions. 

5.2 Volume singularity analysis 

An useful item is to compute the volume time evolution. For that, 
we introduce the time variable ^ defined by 

Adding the three differential equations for a, (5 and 7, and using the 
constraint equation ©, one generally gets the following result for 
u = \n{V) = a + /? + 7: 

u^^ = -2ep + 3C/e^"/2 (9) 



Hi - "4 



-l + (3l+-fl +2F(0) + 3C/e'"/2 = ... (10) 



2u 



al + /?! + 7| + ^e2"C/(0) - ^(^0 + ^) (11) 
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where a subscript ^ means a derivative with respect to ^. This last 
equation can help us to predict cases that are not covered by the 
direct energy analysis. 

In the vacuum Universe, equation Ullj) reduces to 



u 



-2e 



p 



and we have besides the relation: 



(12) 



(13) 



Thus in the vacuum case if the volume starts to decrease, it will 
decrease forever according to equations ()12II13() . indicating a singu- 
larity. 



In presence of the massive scalar field we are back to the full equa- 
tion pij) and we can state the following: 

If [/ < and uo + I > 0, and if the volume starts to decrease, 
it will decrease again forever. The same conclusion applies about its 
possible increase but also for u^(j)^ going to zero with time since we 
have the relation: 

u^^ - n| = V^cp'^dliLnV) + u^c/.^ (14) 

Thus now the quantity u^^ — ti| is not directly connected to the sign 
of the second derivative in time t of Ln{V) due to the coupling term 
in u^(f)^ and to the definition of ^. Hence the conclusions valid in 
vacuum should be adapted including the relative sign of and <j)^. 

• For u^<j)^ > the conclusion remains the same than above 
about a possible volume singularity. 

• For u^(j)^ < the full dynamics should be solved in order to 
compare the relative magnitude of the negative terms in the 
right hand side (rhs) of equation Q and the positive term 

However there is not always a volume singularity. 
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5.3 General predictions 

Our numerical studies have shown that there is not always a sin- 
gularity: in particular the Universe volume can be oscillating with- 
out going to zero. It then exhibites some local minima correspond- 
ing to some bounces Following jl4j . the chaos condition reads 
now h = —H((j)) = P(j,V'^ < while the weak energy condition is 
P<i> + P<f> = —'^Ec{4>)/V'^. Hence it comes: 

• When the volume singularity exists i.e F — > for infinite r 
time, for U bounded in time UV^ and —H — > —Ec{4>) 

(unless the factor ^ in Ec vanishes faster than UV'^). 

Thus there is chaos for < ~f whatever the sign of U, the 
weak energy condition being violated and the energy density 
being negative. There is no chaos for > — |, the weak energy 
condition being respected and p^f) > 0. 

For UV'^ not bounded with time no direct conclusion can be 
drawn. 

• When there is volume singularity or not: 

a) for U > and > — | there is always no chaos, the weak 
energy condition is respected and /j^ > 

b) for U < and cjq < — | there is chaos with violation of the 
weak energy condition and < 

c) for U{ll!o -|- |) < there is no direct conclusion. The weak 
energy condition is violated if u; < —3/2, whatever the sign of 
the potential whereas the sign of p^ is undetermined. 

• Without volume singularity: 

In this case we have a necessary condition that comes from 
u^^ — u| — > which can be split into 2Ec{(j)) — Ep{4>) — 
u^fj)^ > 0. But nothing else can be said. 

5.4 Applications 

Let us make applications for the three specified potentials. 



19 



5.4.1 case 1 : U = 2K 



One concludes from above and also from the constraint equation (jHJ 
that for A > and > — | there should be no chaos in the ap- 
proach to the singularity. While for the case A < and ljq < ~| 
one should observe chaos. This is nicely confirmed by the numerical 
simulations (see figure 

As explained in 111 H 4> oc e'^'^^ , with a constant. Hence 4'(t)~^ = 0o 
and, when there is a singularity (the volume is going there to zero), 
the scalar field kinetic energy will always dominate the potential en- 
ergy on the long time. Thus, in this case, the conclusion for the 
chaos remains the same with respect to the value of loq with respect 
to — |, whatever the sign of A. 

Undetermined situations about chaos occur when there is no singu- 
larity (the volume may be in expansion) and when A(3 + 2ijj < 0). 
Note that in some cases each variable ai may be chaotic while their 
sum (related to the volume logarithm) could be regular. 

5.4.2 case 2: U = (j)^ 

We first consider the case of a real k. Again from the Hamiltonian 
constraint equation (jSJ, one concludes that for u>q > one should 
have no chaos, since the fact to take k real imposes (p > 0. For 
Wo < — I we have no net conclusion since the signs are opposite for 
Ec{(t>) and Ep{(f)). 

Now we consider integer values for k. For k even integer, we expect 
the same behaviour than for k real with (p > 0. But for k odd, we can 
take (j) < and thus a different behaviour of the system is expected 
: this is well confirmed by the numerical simulations (see figure El, 
showing that the notion of energy threshold still holds here. 

5.4.3 case 3: U = e'=* 

Again we are led to the same conclusion than above: there is no 
chaos in the case tA;o > — | whatever the sign of k since the potential 
is positive. However no clear conclusion can be further drawn when 
cjo < — I since 3 -|- 2u>q and U have opposite signs. 

Nothing can be found in the other cases. 
When there is a singularity with the volume going to zero, we try 
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to estimate the value of e^"?7 from equations (|9l lllj) . For instance, 
for A; < and (f) being an increasing function of time, the potential 
term will go to zero and the kinetic energy will dominate at late time 
(idem for the case k > ii (p can decrease with time). Then we get 
the same conclusion as stated above (for a bounded potential energy 
with time). 

6 Discussion 

In this work we studied the dynamical behaviour of the Bianchi type 
IX model in General Relativity with a massive scalar field. We 
consider three forms of potential. A first question was to know if, 
as for the massless Brans-Dicke theory in the Einstein frame, this 
behaviour is ruled only by some control parameters or if the initial 
conditions also play a role. A second question was to know if a scalar 
field always kills the singularity oscillatory approach. 
To answer to the first question, we began to identify three different 
forms of asymptotical behaviours which may be divided in subcases: 

1. The Universe expands in all the directions with > 0. This 
expansion may be singular, i.e. infinite and anisotropic for a 
finite proper time, or lead the forever expanding Universe to 
isotropy. 

2. The Universe contracts monotonically and the curvature terms 
are negligible. Then the solution of the field equations may be 
kasnerian if the potential of the scalar field is also negligible. 
In this case the metric functions behave like power laws of the 
proper time. It may also be non kasnerian if the potential is 
not negligible. In any case, one has > and 3 + 20; > 0. 

3. The Universe oscillates without going to a singularity with 
Pip <0 and 3 + 2a; < 0. 

Contrary to the Brans-Dicke theory in the Einstein frame with a 
massless scalar field, the metric functions behaviour does not de- 
pend only on the control parameters but also on the initial condi- 
tions. However, for each theories, there are some general properties 
which are ruled only by the control parameters. In particular when 
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one considers a theory with a cosmological constant: its asymptotical 
behaviour seems only related to the initial conditions via a special 
value Aq > of this constant when 3 + 2cj > 0. Otherwise, it only 
depends on the signs of A and 3 + 2u;. This weak dependence with 
respect to the initial conditions occurs because oo and A are some 
constants. This theory is thus probably the less initial conditions 
depending one when we consider a massive scalar field. 
For the other forms of potentials which depend on 0, the distribu- 
tion of the metric functions asymptotical behaviour in the control 
parameters plane shares some common properties with the cosmo- 
logical constant one: 

• Oscillatory behaviours only occur when lo < —3/2. Then the 
weak energy condition is violated. The test of chaos resting 
upon the Hamiltonian formalism developed in [T^ and used in 
the previous section goes in this way. 

• Spatially open Universes at late and early times only occur 
when < -3/2 and [/ > 0. 

• Closed Universes at late or early times only occur when a; > 
-3/2. 

• Continuously oscillating Universes only occur for u) < —3/2 
and U <0. 

However, on the contrary to the cosmological constant case, the 
repartition of the expanding or oscillating Universes in the control 
parameters planes restricted by the above conditions strongly de- 
pends on the initial conditions. This is due to the dependence of the 
potential with respect to the scalar field. 

To our knowledge, this is the first time that figures such aslSHHlare 
obtained. They show a new way to study numerically the asymp- 
totical behaviour of a Bianchi type IX model in General Relativity 
with a massive scalar field. It is a new possibility to determine its 
properties in the same spirit as the Brans-Dicke theory in the Ein- 
stein frame with respect to the Brans-Dicke coupling constant uj. 
A similar numerical study was performed in ^2] with the potential 
U = m^cf)^. The sets of initial conditions leading to bouncing solu- 
tions was analysed using a two dimensions slices through the initial 
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conditions space. 

Another common property for the theories of these papers is that 
an oscillating behaviour always needs a negative scalar field energy 
density, does not respect the weak energy conditions and never lead 
to a singularity. An important question is thus to know if an oscil- 
latory singularity approach is always possible despite the presence 
of a massive scalar field. To answer, we rewrote the field equations 
and looked for a particular theory whose scalar field energy density 
becomes sufficiently small to be dynamically negligible. It is what 
happens in General Relativity with a perfect fluid for which the sin- 
gularity approach is oscillatory. We found such a theory and, thanks 
to the potential which allows distinguishing between a violation of 
the weak energy condition and the presence of a negative scalar field 
energy density, we showed that a singularity oscillatory approach 
was possible even with a massive scalar field having a positive en- 
ergy density but always violating the weak energy condition. 
A next step, would be to prove that it also arises when diverges, 
which is the behaviour we are waiting for any kind of matter in the 
neighbourhood of a singularity (if we do not take into account any 
quantum effect). It should be possible since, in the field equations, 
Pff, is multiplied by the comoving 3-volume which should then dis- 
appear. It would also be interesting to test nonminimally coupled 
scalar field or other forms of potential such as the supergravity one. 
Moreover some Painleve tests could be used to analyse the singular- 
ity from the integr ability point of view |24[ I25j . 
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A The Belinskii, Khalatnikov and Lifshitz (BKL) 
approximation 

The BKL approximation ;1D allows understanding easily with heuris- 
tic arguments why the metric oscillates or not. Explanations. 

A.l In General Relativity 

The BKL approximation allows a qualitative description of the met- 
ric functions behaviour for the Bianchi type IX model in the vicinity 
of the singularity. 

One assumes that initially the metric functions are very small. In 
this case, curvature terms appearing in the field equations are neg- 
ligible. These equations can then be approximated by those of the 
Bianchi type / model whose solution is kasnerian: the metric func- 
tions behave as powers of the proper time t such as their exponents pi 
have the properties Y^^=iPi = Z]?=iPi = 1- These equalities imply 
that two of the metric functions decrease while the third one grows. 
It follows that after a finite time, one of the curvature terms, say e", 
is not negligible any more. Then, the BKL approximation is invalid 
and the field equations approximated by that of the Bianchi type // 
model. However, asymptotically this model tends to a Bianchi type 
/ one and thus, the BKL approximation becomes valid again. Re- 
peating this process, the Bianchi type IX model is thus successively 
approximated by some Bianchi type / and // models and the metric 
functions oscillate in the vicinity of the singularity. The interval of 
time between two oscillations is called a Kasner epoch. The time 
during which two same metric functions oscillates is called a Kasner 
era. 

A. 2 In General Relativity with a massless scalar field 

The only difference with the previous case is that from now on, when 
the field equations are approximated by a Bianchi type / model, 
one has YJi=iPi = 1 but J2i=iPi = 1 - + 3/2) where a 
constant. Consequently, when uj > —3/2, all the exponents pi can 
be positive. Then the metric functions are all decreasing and the 
BKL approximation remains valid forever: the oscillations stop. The 
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energy density of the scalar field is positive and it respects the weak 
energy condition. It is the case considered in the paper On 
the contrary, when to < —3/2, the exponents are never all positive. 
Consequently, the BKL approximation is not stable. It is invalid 
after a finite time. In this case, the energy density is negative and 
the weak energy condition is not respected. 
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Figure 4: Asymptotical behaviours of the Bianchi type IX model with a theory 
defined by oj = a)o (in ordinate) and U = 4''° (k in abscissa) when k is real. The left 
graphs represent the early times and the right graphs, the late times. 
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Figure 5: Asymptotical behaviours of the Bianchi type IX model with a theory 
defined by a; = wo (in ordinate) and U = {k in abscissa) when k is an integer. The 
left graphs represent the early times and the right graphs, the late times. 
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Fi gur6 6 1 Asymptotical behaviours of the Bianchi type IX model with a theory 
defined by t<; = wo (in ordinate) and U = e*"^ (fc in abscissa). The left graphs represent 
the early times and the right graphs, the late times. 
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